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Abstract 

We investigate the quantum properties of the weh-known process of sum frequency generation, 
showing that it is potentiahy a very useful source of non-classical states of the electromagnetic field, 
some of which are not possible with the more common techniques. We show that it can produce 
quadrature squeezed light, bright bichromatic entangled states and symmetric and asymmetric 
demonstrations of the Einstein-Podolsky-Rosen paradox. We also show that the semiclassical 
equations totally fail to describe the mean-field dynamics when the cavity is strongly pumped. 

PACS numbers: 42.50.Dv,42.65.Ky,03.65.Ud,03.67.Mn 
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I. INTRODUCTION 



The nonlinear optical process of sum frequency generation (SFG), also known as non- 
degenerate upconversion or frequency summation JlJ, has been known and investigated for 
a number of years, going back at least to theoretical investigations in the seminal paper of 
Armstrong et al. {2! and experimental realisation by Bass et al. jsj. The process has many 
uses outside the quantum optics community, such as, to name only a few, surface vibrational 
spectroscopy of molecules 4|, two-dimensional vibrational spectroscopy [si], studies of liquid 
interfaces 0, 0] , low noise optical tomography 8], and the investigation of powder supported 
catalysts |9|. 

However, despite the wide number of what we may term classical uses for this process, 
very little attention seems to have been paid to its quantum properties, especially in terms 
of entanglement and quadrature squeezing. For the intracavity process, a notable exception 
is the theoretical work of Eschmann and Reid [ic]], who predicted sub-Poissonian photon 
statistics in the high frequency mode and in the sum of the two low frequency modes. More 
attention has been paid to type II second harmonic generation, where the non-degeneracy is 
in polarisation rather than in frequency, for example work by Jack, Collett and Walls 
and by Andersen and Buchhave 12|. Four-wave mixing is also a process which can lead to 
sub-Poissonian fluctuations in intensity sums [13.|], as well as a source of entangled beams 
of either photons [IJ] or atoms This leads us to expect that sum frequency generation 
may also be a source of entangled output beams. 

In this work we will first introduce the full Hamiltonian of the intracavity process before 
examining the properties of the interaction Hamiltonian considered in isolation. We will 
then analyse the outputs of the full intracavity system in terms of squeezing and entangle- 
ment. We will show that SFG is in fact a versatile source of entanglement resources which 
can be easily tuned to entangle beams which have large frequency differences and produce 
entanglement between output beams at different intensities. The tunability of the process 



also allows us to predict t 
and asymmetric steering 



lat the same device could be used to demonstrate both symmetric 



16|, which may have applications in the field of quantum cryptog- 



raphy. We will show that this process is potentially very useful for the field of continuous 
variable quantum information. 
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II. SYSTEM AND HAMILTONIAN 



The basic interaction is that of a photon at ooi combining with a photon at 002 to produce 
a photon at cu^ (= UJ1 + UJ2), mediated by a second order, x^'^\ nonhnearity. (For an accessible 
description of advances in the use of x*'^'' materials, see Hanna [l7|.) The full Hamiltonian 
describing this interaction inside a triply resonant optical cavity and the interaction of the 
cavity fields with the external fields may be written as 

where the interaction Hamiltonian in the appropriate rotating frame is 



(2) 
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Hint = i^fi 

the pumping Hamiltonian is 

"Hpump — ^^^^ [^i^l ) 
1=1 

and the reservoir damping Hamiltonian is 

H,e. = ^E[f4 + fla,]. (4) 

i=l 

In the above, is the bosonic annihilation operator for the mode at frequency Ui, k repre- 
sents the effective second order nonhnearity, the are the classical pumping laser amplitudes 
at the respective frequencies, and the Fj are the annihilation operators for bath quanta, rep- 
resenting losses through the cavity mirror. 



III. EQUATIONS OF MOTION AND HAMILTONIAN DYNAMICS 

We will begin by giving the equations of motion which result from considering the interac- 
tion Hamiltonian in isolation and then add cavity loss and pump terms below, in section IIVI 
as the addition of these is a trivial matter. We note here that the approach we use is not 
expected to give completely accurate predictions for the system operating in the travelling 
wave configuration as it does not account for such physical features as dispersion in the non- 
linear medium. However, it does give an approximate idea and is useful for understanding 
the quantum dynamics which the Hamiltonian makes possible. Following the usual pro- 
cedures jisl . [19I . we may map the interaction Hamiltonian onto a Fokker-Planck equation 
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for the P-representation pseudoprobability distribution of the system [20|, l2l| . Making the 
operator correspondences di ^ ai and d\ ^ a*, we obtain 



dP 



d 



d 



d 



oai oal oa2 



d 



d 



K 

+ 2 



92 



+ 



"3 + 



+ 



92 



5a? 



■Q!]^Q!2 



P(ai,a*,t). 



(5) 



(9a;i(9a2 da2dai ) \da\dal daedal 
We immediately see that the diffusion term of the above Fokker-Planck equation is not 
positive-definite, so we will use the positive-P representation, which at the cost of doubling 
the dimensionality of the phase-space, allows us to map the resulting Fokker-Planck equation 
onto a set of stochastic differential equations 221]. Making the changes a* — ^ af and noting 
that af = a* in a distributional sense, we find the set of Ito calculus (23[ equations, 

dai 
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where the Gaussian noise terms have the correlations 



vAt) = o, vj{t)vk{t') = SjkS{t-t'). 



(6) 



(7) 



We note here that the stochastic positive-P representation equations are an exact mapping 
from the Hamiltonian and allow the calculation of any normally-ordered operator moments 
through an averaging process over trajectories, such that 



(8) 



in the limit of a large number of stochastic trajectories, with the bar representing a classical 
averaging. We also note that this relationship is only valid where there are no divergence 
problems with the stochastic integration, and that there were none in any of the results 
presented here. 



It is worthwhile examining the dynamics predicted by the stochastic integration of Eq. [6] 
as this will provide some insight into the full system where we will consider a doubly-pumped 
triply resonant cavity. We note here that we will not attempt to analytically solve the mean- 
field equations found by dropping the noise terms in Eq. [HI as this procedure is of limited use 



in second harmonic generation 



24 



25 



26l |. which differs in having the two lower frequency 



modes as identical. Hence we consider that linearisation of these equations is unlikely to 
be accurate after some shortish interaction time for either the mean fields or the quantum 
correlations. We will use a linearisation procedure in section IIV| where this is much more 
useful as long as we are careful with the regime of viability. 
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FIG. 1: (colour online) Mean intensities from an average over 2.67 x 10^ stochastic trajectories of 
the positive-P equations ([6]). The horizontal axis is a scaled dimensionless time, C = «^|ai(0)|t. All 
quantities plotted in this and subsequent graphs are dimensionless. 

The first quantities we calculate are the mean intensities as a function of interaction 
time, as shown in Fig. [T]for k = 0.01 and the initial conditions ai{0) = 02(0) = lOOO/y^ 
and 0:3(0) = 0. Both the pumping fields are treated as coherent states, which is reasonable 
for stabilised lasers operating above threshold. We see that the dynamics are reminiscent 
of those obtained from a similar treatment of second harmonic generation, showing an al- 
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most complete conversion to the higher frequency, followed by partial reconversion to the 



two lower frequencies 



25| . We note here that a previous semiclassical analysis for equal 



intensities in the two low frequency inputs showed complete conversion and no subsequent 
reconversion 27|, processes which can be explained as due to quantum fluctuations and the 
discrete nature of the electromagnetic field [25|, l28j. Stochastic integration using the phase- 



space representations automatically includes these effects and allows for different quantum 



states of the input fields to be simulated, such as thermal 



29| and number states [30|. 
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FIG. 2: (colour online) The X quadrature variance for the field at and the Fano factor for 
the intensity sum of the two low frequency fields. A value of less than one signifies squeezing or 
sub-Poissonian statistics respectively. The two low frequency fields show excessive noise in both X 
and Y quadratures and exhibit super-Poissonian intensity fluctuations. 

As our motivation for investigating this system is its utility for producing entangled 
and other quantum states of electromagnetic fields, we will now look at these quantum 
properties. The first of these is single-mode squeezing, for which we first need to define 
quadrature operators. We will define a general quadrature of the field at a given quadrature 
angle as jl8l | 

X^{e) = aje-'^ + a]e'^, (9) 
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and use the shorthand Xj{0) = Xj and ^j(f) = Yj. With these definitions, a squeezed state 
is one for which V{X{9)) < 1. In Fig. [2] we show the results of stochastic calculations for 
the variances of both the high frequency X^ quadrature and the sum of the low-frequency 
intensities. For the latter, which has been predicted to give sub-Poissonian statistics in the 
intracavity case pLOj], we have plotted the Fano factor, defined for this 

V{Ni + N2) 



F{Ni + N2 



(10) 



N1 + N2 ' 

so that any value below one means that the combined mode has less intensity fiuctuations 
than a coherent state. Comparing the two curves in Fig. [21 we see that the squeezing in X3 
essentially continues until the upconversion turns into downconversion, whereas the noise 
suppression in the intensity sum disappears as the low frequency intensities approach their 
minima. The other quadratures, the single-mode intensities and the intensity differences 
all exhibit excess noise as the interaction proceeds. This ability to produce quadrature 
squeezed light at the sum frequency could be useful in, for example, super resolution optical 



measurements 



3l|. 



The correlation between the intensities of the two low frequency fields suggests that 
they may exhibit bipartite entanglement. We examine this using the criteria for continuous 



variables developed by Duan et al. 



32l | and Simon 33|]. In the present case where the two 



pump modes have equal intensities, it is sufficient to violate one of the inequalities, 

V{X^ ± X2) + r (fi T %) < 4, 



(11) 



to establish that the two low frequency modes are entangled. It is known that entan- 

n 

gled states are a superset of states which can be used for steering [16!] and that this can 
be demonstrated in the continuous variable case by violation of the Reid EPR (Einstein- 
Podolsky- Rosen) inequalities SJ, l35|] . In this case the inequality can be written as 



V,rrfiXj)Vir.fiYj) > 1, 



(12) 



where 



Vinf{Xj] 
VinfiYj] 



V{X,) - 



V{Xj,Xk 
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V(Y^,Yk] 



V{Yk) 
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FIG. 3: (colour online) The Duan-Simon correlation, +X2) + V{Yi — Y2), and the Reid EPR 

correlation, Vinf{Xj)Vinf{Yj), for the two low frequency fields. With symmetric inputs, the EPR 
correlation is equal for j = 1 and 2. Note that the Duan-Simon value has been divided by four so 
that both entanglement and the EPR paradox are demonstrated for values less than 1. 

with violation of the inequality showing that modes j and k are entangled in the EPR 
sense. As can be seen in Fig. |3l where we have divided the correlation for the sum of the 
X quadratures and the difference of the Y quadratures by 4, these two measures give an 
unambiguous demonstration that the two modes become entangled by the interaction and 
are in fact entangled in the strong sense required by the concept of steering. Having now 
established the properties of the Hamiltonian dynamics, we will turn our attention to a more 
quantitative analysis of the intracavity case. 



IV. INTRACAVITY DYNAMICS 



In the intracavity case we are interested in correlations of the steady-state fields, as these 
are what are normally measured using homodyne techniques. To analyse these theoretically, 
we must add pumping and loss terms to the equations given above ([6]). Making the usual zero 
temperature Born and Markov approximations 18|] for the interactions with the reservoirs. 
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we find 
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73a^ - Kafat, (14) 

where the jj are the cavity loss rates at the respective frequencies, the ej are coherent 
pumping terms, and the noise terms have the same correlations as in Eq. El 

In the intracavity case it is often possible to find the appropriate fluctuation spectra via 
a process of linearising the fluctuations in the variables about their steady-state, classical 



solutions 



This enables us to write equations for the fluctuations as a multivariate 
Ornstein-Uhlenbeck process, from which it is particularly simple to extract the appropriate 
spectral correlations. There are however, two caveats which must be considered here. The 
first is that the steady-state solutions of the classical equations (Eq. [THwith the noise terms 
dropped) must be the actual mean-field solutions, and the second is that the drift matrix in 
the resulting equation for the Ornstein-Uhlenbeck process must not have any negative real 
parts to its eigenvalues. These conditions are known to occur with intracavity parametric 



processes, with 



eration 



36 



37 



;he first being violated in the self-pulsing regime of second harmonic gen- 



38| and the second being a problem at the critical operating point of the 



optical parametric oscillator 39| . In order to check the stability we first turn to solving the 
steady-state equations for the mean fields. 

Although in general not possible, we can find analytical solutions for the steady states 
in the case where we set the pumping strengths equal, with ei = €2 = e, and the two low 
frequency loss rates also equal, with 71 = 72 = 7- Although achieving equal loss rates 
may not be so simple in the laboratory, it is useful here in order to gain some insight. 
With these values set equal, the two low frequency modes will have equal amplitudes with 
a^* = = a***, which leads to a cubic equation for Og** (note that we will now drop the 
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superscripts for notational convenience) 



73K^a3 - 2773Ka3 + 737^03 + ft^e^ = 0, 



(15) 



along with a simple expression for a, 

a = . (16) 

7 - 

We note that there is no divergence here in the value of a since the high frequency field 
amplitude is negative in the steady state. Without setting these loss and pumping rates 
equal, it is much more difficult to find analytical solutions, although we can proceed via 
numerics, which we will do later in order to examine the effects of asymmetry. In the 
present case where the cavity is triply resonant, 03 will be the real solution of Eq. [151 This 
is found as 

(17) 



as = — 
^ 6 



47 16V3^^73 
K ^ 73fi:2 



where 



,2,^5^2 



.nl/3 



(18) 



We now decompose the variables into their steady-state classical values and fiuctuations 
around these, 



a 4 



(19) 



and find the equations of motion for the fiuctuation vector. 
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6ai, 6ai ^ 6ai^ ^(yti ^^^i 



(20) 



as 



5X = -Adt + BdW, 



(21) 



where A is the drift matrix (remembering that the aj are now to be read as the steady-state 
values) , 
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10 



dW is a vector of real Wiener increments and B is the diffusion matrix, 



B 









ij^ 





/Ka£ ^ 
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-i^'^ 























(23) 



In parameter regimes where the matrix A has no negative real part to any of its eigenvalues, 
we may simply find the intracavity spectra via the relation 



S{uo) = {A + iujl)~^ BB'^ i^A^ -iuol) ^ 



(24) 



after which we use the standard input-output relations 40| to relate these to the experimen- 
tally measurable quantities outside the cavity. 
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FIG. 4: (colour online) The regions of stability and instability as a function of e and 73/7 for 
K = 10~^ and 7 = 71 = 72 = 1. In the region to the right of the line there is always at least 
one eignevalue with a negative real part and the linearisation process for calculating spectra is not 
valid. 
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Again setting ei = 62 = e, and the two low frequency loss rates also equal, with 71 = 72 = 
7, so that ai = 0^2 = a and all the intracavity fields are real, we find analytical solutions for 
the eigenvalues as 

Ai,2 = 7 ± /taa, 
1 



A3,4 - 2 

X _ 1 

A5,6 - 2 



7 + 73 + «:a3 ± \J [l - l^Y + 2Ka3(7 - 73) + /«^(a3 - 8"^) 
7 + 73 - Kftg ± v^(7 - 73)^ - 2Ka3(7 - 73) + /t^(«3 - 



(25) 



Since 03 is negative in the steady-state, there is an obvious instability when < — 7/fi:, 
found from the first pair of eigenvalues. Unlike the case of the optical parametric oscillator, 
where a simple expression can be found for the critical pumping value, we find that it is 
easiest here to express this in terms of the low frequency field amplitudes, defining a critical 
field, 

a. = ^. (26) 

We note here that this simple expression will not hold unless we set the pumping and loss 
rates symmetrically, as we did for our analytic eigenvalue analysis. Fig. H] shows how the 
parameter space is divided into stable and unstable regions in this symmetric case, as we 
vary the pumping and the ratio between the high and low frequency loss rates. All the 
spectral results we give will be from the stable region. We will return to the region above etc 
below, using stochastic integration to solve the full equations of motion. In the meantime 
we will examine the quantum correlations in the stable region of the parameter space. 

V. ENTANGLEMENT PROPERTIES 

In order to examine the entanglement and squeezing properties of this system, we must 



go beyond the intensity fiuctuations calculated previously [10|] and calculate the phase de- 
pendent quadrature correlations given above in Eq. [11] and Eq. [121 The first of these which 
we show, in Fig. [5], is the output spectral variance of the high frequency amplitude quadra- 
ture, X3. We see that, as in the travelling wave case, this is squeezed and that the degree 
of squeezing increases as we approach ac- As shown in Fig. [HI this system is also a good 
source of bichromatic bipartite entanglement between the two low frequency fields. We see 
that the degree of violation of the Duan-Simon inequality also increases as we approach 
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although not as markedly as does the high frequency squeezing. Consistent with the 
fact that SFG suppresses the noise in the sum of the two low frequency intensities, unlike 
parametric downconversion which suppresses the noise in the difference, we find that the 
entanglement is signalled by the correlation V{Xi + X2) + — ¥2)- We did not find any 
evidence of tripartite entanglement in this system, at least for the parameter ranges we have 
investigated. 
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FIG. 5: (colour online) The spectral variance of the high frequency quadrature, X3, for k = 10~^, 
71 = 72 = 7 = Ij and 73 = 10. The solid line is for ei = £2 = 200, giving ai = 02 = 0.25ac, 
the dash-dotted line is for ei = 62 = 400, giving ai = 02 = 0.61ac) and the dashed line is for 
ei = €2 = 600, giving ai = 02 = 0.95ac- 



A. Symmetric and asymmetric bipartite steering 



In a recent article 16|, Wiseman et al. discussed the concept of steering, which had 



originally been introduced by Schrodinger 4l|, |42] in the context of the EPR paradox 43|. 



Wiseman et al. explain this concept in terms of a bipartite state prepared by Alice, who 
then sends one part to Bob, with this process being repeatable. After measuring their own 
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FIG. 6: (colour online) The Duan-Simon entanglement correlations for the low frequency modes 1 
and 2 and the same parameters as Fig. [5l The two quadratures which violate the inequality (less 
than 4) and hence signify entanglement are X1+X2 and Y1 — Y2. The solid line is for ei = e2 = 200, 
giving ai = a2 = 0.25ac) the dash-dotted line is for ei = £2 = 400, giving ai = 02 = 0.61ac, and 
the dashed line is for ei = €2 = 600, giving ai = a2 = 0.95ac. 

parts, they communicate classically, with Alice trying to convince Bob that the prepared 
state is entangled. If Bob cannot explain the correlations using any local hidden state (LHS) 
model, the state must be entangled. Schrodinger used the concept of LHS to say that Bob's 
system could have a definite state before measurement, even though this actual state would 
be unknown to Bob. He introduced the concept of steering to describe how Alice could 
affect Bob's state via her choice of measurement basis, but expected that this would never 
be seen experimentally. 

For the purposes of this paper j_jt is interesting to note that a demonstration of the EPR 
paradox using the Reid criteria |34|, as was first done by Ou et al. 4J] using parametric 
downconversion, is equivalent to a demonstration of steering. In Fig. [7| we show that this 
system does allow for demonstrations of steering, with a degree of violation of the Reid 
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FIG. 7: (colour online) The spectral EPR product of the inferred variances for modes 1 and 2, 
and the same parameters as Fig. [3 The solid line is for ei = £2 = 200, giving ai = 02 = 0.25ac, 
the dash-dotted line is for ei = €2 = 400, giving ai = 02 = 0.61ac, and the dashed line is for 
ei = €2 = 600, giving ai = 02 = 0.95ac. 

inequality f|T2l) that increases as the pumping is increased so that etc is approached. For the 
symmetric inputs and parameters we have chosen here, the steering is totally symmetric, 
with either Alice (ui) being able to steer Bob ((^2) or Bob being able to steer Alice. However, 
in the work of Wiseman et ai, the reader is left with an open question as to whether there 
exist asymmetric states that are steerable by Alice but not by Bob 16|. It seems intuitively 
obvious that these type of states would not be produced in any normal downconversion 
processes, which are inherently symmetric in their production of correlated pairs of photons. 
The freedom we have with the present system, where the pumping rates and mirror losses 
can be different at the two low frequencies, leads us to expect that it should be a good 
candidate for the exhibition of asymmetric steering. Indeed, as shown in Fig. [HI it is in 
principle possible to arrange the cavity loss rates and pumpings so that this is seen. In the 
figure, we have used the labels EPRij to represent the product Vinf{Xi)Vinf(Yi) and have 
chosen a regime where the linearisation process is valid. We also note here that previous 
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coupled systems which have been predicted to produce entangled outputs via evanescently 



coupled intracavity x^^^ and x^^^ processes could also be arranged asymmetrical 
may also be good candidates for the production of asymmetric steering 3^, 45 . 



y and thus 



46 



47 



48|. 



Having shown that asymmetric steering is possible, we may pose another open question as 
to what it may prove useful for, with one way quantum cryptography being one suggestion, 
although further development of this theme is outside the scope of this paper. 
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FIG. 8: (colour online) Demonstration of asymmetric steering for the parameters 7i = 1, 72 = 
40, 73 = 2, ei = 400, £2 = 2400 and k = 0.01. The correlation EPR12 shows that 2 can steer 1, 
while EPR21 shows that 1 cannot steer 2. 



B. Quantum dynamics 

In the unstable parameter regime as shown in Fig. H] we must resort to stochastic in- 
tegration of the full equations dill), without linearising about the semiclassical solutions. 
What we find is qualitatively different from the behaviour in the stable regime, as can be 
seen in Fig. |9l For pumping strengths for which the semiclassical solution has a < ac, the 
full quantum solutions converge to the semiclassical predictions, as they do in the transient 
regime of Fig. [91 However, when the cavity is pumped more strongly, the fields change 
radically after a short time, with the low frequency fields increasing and the high frequency 

16 



intensity decreasing as the process of downconversion back into the two low frequencies 
becomes dominant. While interesting from the point of view of being a system where the 
semiclassical solutions do not capture the dynamics of the mean fields, we found no evidence 
of steady-state entanglement or quadrature squeezing in our investigations of this regime. 
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FIG. 9: (colour online) Intracavity intensities predicted by averaging 2.2 x 10^ stochastic trajectories 
of the full positive-P representation equations, for 71 = 72 = 1, 73 = 10, ei = £2 = 10'^ and k = 0.01. 
The straight lines are the semiclassical steady-state predictions for the same parameters. 



VI. CONCLUSION 

In conclusion we have shown that the system of intracavity sum frequency generation is 
potentially a versatile and tuneable source of squeezed light, bright bichromatic continuous 
variable bipartite entanglement, and two mode Einstein-Podolsky- Rosen states. With recent 
advances in the fabrication of nonlinear crystals and with the techniques necessary for what 
we may term its classical uses very well established, this process may become a powerful 
tool in the optical arsenal of quantum information. The potential ability to run the system 

17 



in a highly asymmetric manner means that it could be used as an instrument for research 
into the fundamentals of quantum mechanics in ways that are not possible with the familiar 
optical parametric oscillator. The ability to choose and combine a range of frequencies, with 
the combination frequency being quadrature squeezed, may prove useful in super resolution 
measurements. 
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